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Abstract 

Quantum fields exhibit non-trivial behaviours in curved space-times, especially around black 
holes or when a cosmological constant is added to the field equations. A new scheme, based 
on the Wcntzcl-Kramers-Brillouin (WKB) approximation is presented. The main advantage of 
this method is to allow for a better physical understanding of previously known results and to 
give good orders of magnitude in situations where no other approaches are currently developed. 
Greybody factors for evaporating black holes are rederived in this framework and the energy 
levels of scalar fields in the Anti-dc Sitter (AdS) spacetimc arc accurately obtained. Stationary 
solutions in the Schwarzschild-Anti-de Sitter (SAdS) background are investigated. Some im- 
provements and the basics of a line of thought for more complex situations are suggested. 
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1 Introduction 



The study of some of the most intriguing phenomena in our Universe, such as the evaporation 
of microscopic black holes or the parametric amplification of primordial quantum fluctuations, 
is performed within the framework of quantum field theory in curved space-time. Although 
well-established and reliable, as it does not require any quantization of the gravitational field 
itself, this theoretical framework leads to highly non-trivial partial differential equations whose 
resolution is, in most cases, based on numerical investigations. The Wentzel-Kramers-Brillouin 
(WKB) approximation and its extension to semi-classical methods (introduced by Kramers m. 
Young & Uhlenbeck j21 and rcaddressed by Langer can be extremely usefuU to describe 
those physical systems. 

The WKB method is used to approximately solve the following equation : 

{■§2+p'(o)m^o (1) 

where the real function can be either positive or negative valued. Nevertheless, WKB wave 
functions ip are exact solutions of another differential equation : 

which restricts the validity of the approximation to regions where 

\P^ \ ^ 

It means that if the momentum is a linear function of ^, the WKB wave function becomes an 
exact solution of Eq. . The WKB approximation can therefore provide very accurate results 
if the considered momentum squared is, depending on the potential, dominated by the constant, 
linear and quadratic terms in its Taylor expansion. 

Recently, this approximation has been used to derive many important results about quantum 
fields propagating in a curved background. The Hawking evaporation process 0j, as well as 
the inflationary cosmological perturbations [SjE], have been investigated within this scheme. 
In Ref. jj, the Hawking temperature for a Schwarzschild black hole was rederived whereas, in 
Ref. the spectral index and the amplitude of the inflationary perturbations were recovered 
nearly exactly. Finally, the WKB method has also been extensively used in order to determine 
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the quasi- normal frequencies of black holes (see Ref. [7] and Ref. ^J)- The WKB tunnel 
probability used in the computation of the quasi-normal frequencies has also been applied to 
estimate the greybody factors for gravitons emitted by a D— dimensional Schwarzschild black 
holes [ini. 

In this article, we extend the idea of the WKB and semi-classical methods to the propaga- 
tion of massless, minimally coupled, scalar fields in a static and spherically symmetric curved 
background, called Schwarzschild-like hereafter. Section 2 is devoted to the derivation of the 
radial part of the Klein-Gordon (KG) equation and to its rewriting in a Schrodinger-like form, 
such as Eq. 1^. It is shown that this new equation is the quantum version of a slightly- modified 
relativistic Hamilton- Jacobi equation. The semi-classical scheme is applied to solve the prob- 
lem and to derive the quantization rule and WKB tunnel probability. In Section 3, the discrete 
energy spectrum for massless scalar fields embedded in an Anti-de Sitter (AdS) Universe is 
derived and compared with the exact spectrum In Section 4, the problem of the deter- 

mination of the greybody factors, i.e. the coupling between quantum fields and black holes, 
is investigated. The emission/absorption cross-sections are thus determined for scalar fields 
evolving in a 4— dimensional Schwarzschild background. The resulting WKB cross-sections are 
also compared with the exact ones (and to those obtained in Ref. jni)- In Section 5, the 
propagation of scalar fields in a Schwarzschild- Anti-de Sitter (SAdS) space-time is addressed : 
it is shown that resonances, whose energies and bandwidths are obtained, arise in the state 
density. Finally, in Section 6, we draw up some conclusions and comments about the usefulness 
and reliability of the semi-classical approximation and underline some perspectives to improve 
this method. Some general arguments about its significance to study quantum fields in curved 
space-time are briefly given. 



2 Scalar fields in Schwarzschild-like space-time 
2.1 Equation of motion 

A static and spherically symmetric space-time is described by the general Schwarzschild-like 
metric, totally determined by one function h{r) : 

ds^ = h{r)dt^ - -S^r - r^dn^ (4) 
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where dfl is the sohd-anglc. The propagation of masslcss and minimally coupled scalar particles 
in this backgromid is given by the solution of the KG equation = 0. The field can be 
decomposed into eigenmodes of normal frequency uj and angular momentum i, <f>(t, r, f2) = 
e~'^^'^Y™{n)R{r), the dynamics being totally given by the radial hmction R(r) which satisfies 
the radial part of the KG equation : 



d / , n dR 



2^2 



where £{£ + 1) corresponds to the eigenvalues of the angular part of the KG equation. 

The resolution of this equation strongly depends on the shape of the metric function h 
and it is usually not possible to obtain exact solutions, even in a simple background such as 
the Schwarzschild space-time. However, Eq. ^ can be modified by the use of the tortoise 
coordinate r*, defined as dr* = dr/h{r), and a new radial function U{r) = rR{r). It becomes 
of the Schrodinger type : 

V ' 

which means that the square of the momentum operator id/dr* is equal to — M{r*). We 
interpret Eq. © as the operational version of the classical Hamilton-Jacobi equation = 
p2 _^ ^2^^* -J ^j^gj-g function is identified to the square of a potential-like contribution. 

Except by the author of Ref. |E] j the M function is usually associated with V and not with 
(see, e.g., Ref. Although this does not change the mathematical structure and the 

solutions, we prefer to identify M function with the square of the potential in order to obtain 
a homogeneous classical Hamilton-Jacobi equation in a physically-consistant way. Eq. (jSJ can 
be regarded as the quantum version of the classical equation of motion for a test particle in a 
Schwarzschild background given by : 



\dT J 



r 



where the second term of the left hand side of the above equation is seen as the square of the 
potential. Consequently, the propagation of scalar fields in a Schwarzschild-like background is 



3 



based on a classical problem described by the action S solution of 

It should be noticed that the above equation slightly differs from the relativistic Hamilton- 
Jacobi equation for massless particles evolving in a spatial potential. Particles propagating 
in a spatial potential are described by the Lagrangian C = — 7ti\/1 — — V{r*) whereas the 
considered problem can be seen as a free particle with a mass which depends upon the position 
C = —V{r*)\/1 — v'^. Because Eq. Q has exactly the same mathematical form as Eq. (Q, the 
semi-classical approximation can be easily applied to the propagation of massless scalar fields 
in a Schwarzschild-like background, keeping in mind that the underlying classical problem is 
described by the above-mentioned Lagrangian. 

2.2 Semi-classical approach 

The quantum version of the Hamilton- Jacobi relation is given by the following KG-like equation : 

(-^^^^ Htt) = -f^'^i^ty (9) 

In a Schwarzschild-like background, this equation corresponds to the time and radial part of 
the KG equation, the angular part being separated using spherical harmonics, with ^ standing 
for the tortoise radial coordinate and t for the Schwarzschild time. The U radial function of eq. 
|S1 corresponds to the spatial part of the ip wavefunction in the above equation. For an initial 
event 2 = and a final one J-' — (i/,^/), the ansatz corresponds to the sum over all the 

classical paths given by 

V'(.F,X)=^F„(.F,X)e'^^, (10) 

n 

where Sn{^,2) is the classical action of the n*^ path between T and JF. Performing the power 
expansion up to the linear order leads to : 
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where the n indices have been dropped for simphcity. The resolution is not directly performed 
and an alternative set of equations should be used. The first order relation leads to : 



Considering Eq. H12|) as an equation of continuity, we can define the quadri-current at the 
semi-classical order as \F\^ Pf^, so that : 

-|F|'p^=i(V'9^V^-^^a^^). (13) 

For a stationary problem, i.e d'^ip/dt^ = ^Lu'^tp, it is more convenient to consider a path 
with fixed energy and the Green function of the propagator is preferred. Using the stationary 
phase approximation, it is shown that, as in non-relativistic quantum mechanics, the Green 
function takes the form 



G(^/,C.,-)=5:./^e^4^-(0^^ (14) 
„ V P'^^^f) 

where p„ = corresponds to the classical momentum ^uP- — ^^(0) ^^'^ energy oj and 
does not depend upon the considered path. Althought most non-relativistic quantum results 
can be easily adapted, we have carefully checked, thanks to an explicit calculation of the density 
of states, that the Bohr-Sommerfeld quantization rule in a well-potential is recovered : 

2 /^^^ = (2n+l)7r, (15) 

where n is an integer and ^± the two turning points. Following the techniques of Ref. |il6| (the 
classical action is here allowed to become complex) and using the matching procedure at the 
turning point described in appendix IXI the transmission coefficient at the WKB order is found 
to be 

|.4|'=exp(-2 r%'(C,c.)de) (16) 



where the interval [^_,^-|-] corresponds to the region where energy of the particle is lower than 
the potential. 
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3 Energy spectrum in AdS Universe 



3.1 Exact results 

It was shown in that massless scalar fields, minimally coupled to gravity, propagating in an 
AdS space-time exhibit a discrete energy spectrum. We will briefly rederive those exact results 
using another -less rigourous- approach to compare them with semi-classical computations. For 
scalar particles propagating in a pure AdS space-time, the following Schwarzschild-like metric 
has to be considered : 

where R = (—A/3) is the curvature radius of the AdS space-time. The KG equation is solved 
[S] using a more convenient coordinate system |17[ [T5| : 



r = ^ ; p = 2arctan ( ^ + ^1 - ) - | ; = ^ {dr^ - dp^ - sin^ ip)dn^) . (18) 



With the ansatz (f>(r, p, fi) = e '"''^Yp {Q)R{p), the following radial equation is obtained : 

cos2 (p) d fcos^ ip) dR\ , f^^ _ e{e + i) \ ^ ^ Q ^^g^ 



sin (p) dp \ sin (p) dp J \ sin (p) 
which leads -see appendix E]- to : 

j(<?+l/2,3/2) 



Rip) = sin^ ip) cos^ {p)Pir/J,;i,' (cos {2p)) . (20) 



This solution, based on Jacobi polynomial, leads to the energy spectrum 



'2n + £ + 3\ 



where a zero point energy clearly appears. The curvature of space-time, 1/ R, provides a natural 
unit for scalar particles energy and this quantization can be viewed as periodic conditions since 
the geodesies are closed in the AdS Universe. 

3.2 The Bohr-Sommerfeld quantization rule 

The above results are a direct consequence of the quantum nature of the fields combined with 
the intrinsic curvature of the AdS Universe. Rewriting the radial part of the KG equation in 
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Figure 1: Potential seen by a scalar field in AdS space-time in units of \/\A\ 



a Schrodinger-like form, as mentioned in section 12.11 it is clear that the potential seen by a 
scalar particle is a well-potential for all the values of the orbital momentum but zero. This is 
illustrated in Fig. ^and given by 
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i?2 cos2 {r*/R) Vtan2 {r*/R) 



1) 



+ 2 



(22) 



where the tortoise coordinate r* = i?arctan (r/R) lies in the range [0, 7ri?/2[. As a consequence, 
scalar fields clearly have to exhibit a discrete energy spectrum at variance to what happens in 
an asymptotically flat Universe. Moreover, the Bohr-Sommerfeld quantization method can be 
applied in order to determine the energy levels and compared with the exact results. 



3.2.1 The 



case 



The £ = case is first considered as some analytical calculations can be performed. The two 
turning points are r*_ = and = i?arctan ^2^ ^ ^^'^ potential takes the following 

V2 



form Vo{r*) 



5(?'*), where a Dirac distribution was included so as to take into 



account that the negative r* region is forbidden. When one of the walls is infinitely steep, 
as is the case here, the Bohr-Sommcrfeld quantization rule is slightly modified because of the 
new boundary conditions for the wave function |15| . The quantization rule reads in this case 
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Wq = (2n + I) TT. To account for the distribution involved, a more rigorous definition of this 
quantity shall be used : 



H^o = 21im/ Vtj2 _ v^{r*)dr*. (23) 



We finally obtain the following value for the classical action Wq — tt (loR — V2) which yields 
to the following spectrum for the £ = case : 

(^«,o)„eN = ■ (24) 

The spectrum dependence with respect to the integer n is well recovered and, as a consequence, 
so are the energy gaps. Nevertheless, the zero point energy slightly differs from the accurate 
one : the semi-classical approach provides tJo.o = (3/2 + \/2)/i?~ 2.91/i? instead of 3/i?. This 
shift is due to the specific behaviour of the potential at r* = leading to a high departure 
from any quadratic approximation. It is important to mention that if an infinite wall was not 
added, the quantization rule would follow Eq. (|15f) . yielding to an energy spectrum which mis- 
estimates the zero point energy by 20% because the potential, if viewed as a function of the 
two variables {r*,£), is not continuous at (0,0). 

3.2.2 The i case 

For £ > 0, the analytical computation of Wo{uj) is not possible and numerical investigations 
have to be performed. To evaluate the energy levels, the equation 

Wo{uj,£)-{2n + l)TT = (25) 

is numerically solved. Wq has to be rewritten in a more convenient way using the y ~ tan {r* / R) 
coordinate : 



where y± are the two turning points, determined in appendix lUl 

y± = ^\luj^R^ -2 -£{£+ I) ±^J{2 + £{£+!)- uj^R^f - 8£{£ + 1). (27) 

The numerical results are presented on Fig.[21and Fig.|31where the energy levels u;„^i and the 
relative errors Ato / uj are displayed respectively as a function of the integer n and as a function 
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Figure 2: left 'panel : Energy levels as a function of the integer n. Right panel : Energy levels 
as a function of the orbital quantum number i. All the energies are expressed in units of R~^. 



of the orbital quantum number £. Although slightly underestimated, the energy levels are well 
recovered and the precision of this semi-classical approach increases with n and i : the total 
relative error remains below 5%. For a vanishing orbital quantum number, the result of the 
previous analytical calculation is recovered, as it should be. When £ ^ 0, the energy levels as a 
function of the angular momentum are well described by a linear function with a ^—intercept 
equal to (2n + 3)/i?. Recovering the correct zero-point energy by taking the £ limit lights 
up the particular behaviour of the potential for a vanishing angular momentum and confirms 
the prescription used in the previous section. 

There are two possible sources of error in those numerical computations. The first one is 
due to numerical approximations and was checked to remain totally negligible. The second 
one is related to the Bohr-Sommerfeld quantization scheme which gives exact results up to 
the second order in the potential only. This error decreases both with n and with the orbital 
quantum number i : for a fixed value of n, not only does the energy levels increase with i but 
also docs the minimum of the potential Vmin = (^V^ + ■\/£(£ + 1)^ R^^. This result is expected 
as particles with higher energies are closer to "classicality" . 
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Figure 3: left panel : Total relative error ^j- as a function of the integer n. Right panel : Total 
relative error as a function of the orbital quantum number i. 



4 Greybody factors for Schwarzschild black holes 

For masslcss scalars, the flux emitted by an evaporating black hole can be written as |21| : 

dPN _ 47rcrg([j)a;2 



dudt 



(28) 



The non trivial part mostly lies in the determination of the greybody factors (Jg{^) (related 
with the probability for an emitted particle to escape form the gravitational potential) whose 
computation involves the resolution of the KG equation, given by Eq. (jSJ, in the Schwarzschild 
background. The tunnel probability \A(\^ |22 yields to the greybody factors via the optical 
theorem : 

-.H = E^^^l^^l'- (29) 

Unlike the AdS case, the intricate shape of the metric function describing Schwarzschild 
black holes does not allow for an analytical resolution of Eq. As a consequence, all the 

evaluations of the greybody factors, from usual Schwarzschild black holes to Z?-dimensional 
Gauss-Bonnet ones |22ElEniEnii Eire based on numerical investigations; thought some an- 
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Figure 4: Potential seen by a scalar field propagating in Schwarzschild space-time 

alytical methods can be applied in the high and low energy regimes j27j . However, using the 
WKB method, some approximated results can be obtained without big numerical codes. The 
aim of this study is to derive the greybody factors for 4-dimensional Schwarzschild black holes 
at the WKB order and to compare them with the exact results. 



4.1 WKB wave function and tunnel probability 

For a black hole with a horizon radius rn, we consider the metric : 



r / 



1_ £« 



- r-'dfl'' 



(30) 



In the tortoise coordinates system (t, r*), particles are described by a wave function U evolving 
in the potential V (see Fig. 0J : 



(31) 



As it can be seen from Fig. 01 this potential tends to zero when r* ±oo, making the use 
of the optical theorem valid. We define r^^^ as the tortoise coordinate corresponding to the 
maximum value Vinax- 
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To determine the wave function, the "no- incoming wave at spatial infinity" hypothesis is 
used as there are no sources of particles far from the black hole. For particles with energies 
always greater than the potential, the wave function reads U{r*) — A^^t cxp (^i JJ* p{x)dx^ 
with p = \/ up- — . For particles with energies lower than Vmaxi the space is divided into 
three regions A., B and C, corresponding respectively to ] — oo, rl], [rl , r^] and [r^, +00 [, with 
(rl,7'*) representing the two turning points. With the above boundary conditions, the wave 
fimction admits the following asymptotic regime : U{r*) ~ Ain exp (—icur*) + Aout exp (icor*) 
when r* —00 and C/(?'*) — -Bout exp [iujr*) when r* — > +00. The WKB wave function must 
therefore be read as : 



1 / A — * p(x)dx . i f^^, p(x)dx\ . . . 

A„e ^''1 +Aoute ^"-i ' m region^ 

^c^a jt^ m region O (ozj 



^o^t region C 



with p and p' defined as in section [2. 21 The above solutions show that particles in region A are 
described by both incoming waves, i.e particles emitted by the black hole, and outgoing ones, 
i.e particles reflected on the potential barrier, whereas they are only described by outgoing 
waves in region C. In region B, particles are only described by evanescent waves as there are 
no sources of particles at spatial infinity. Using the procedure described in appendix to 
match the different solutions at the two turning points, it is found that Am = e^^^^^Aout and 
Bout = Aout with r = J^J^ p' (x)dx. This approach leads to a new understanding of the 
evaporation process. 

Unlike the AdS case, it is not possible to inverse the relation r* = f{r) and an explicit 

expression of the potential as a function of r* cannot be found. However, the determination of 

the two turning points as well as r*^^^^ can be performed using r as a variable instead of r* . 

The tunnel probability is given by the ratio of the outgoing flux Tjj ~ i{Udr*U'^ — Wdr*U) at 

spatial infinity to the one at the event horizon of the black hole. Once recxpressed using the 

_ 2 

usual radial coordinate, the tunnel probability Ag is given by : 



A, 



2 1 if Lu > 



max 



„ ^^^^ 

M.) if ,u<Vn, 



^ max 



where r± are the two turning points and Vmax is obtained by computing the value of the 
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potential at Tmax so that 

|rff for ^ = 



,x = < 



(34) 



+ v/9 + 9£2(£ + 1)2 + 14^(^+1) y 



for ^ 7^ 0. 



Although the maximum value of the potential is analytically determined, the two turning 
points as well as the integral involved in the calculation of the tunnel probability have to be 
numerically computed. The above tunnel probability has to be linked with the transmission 
probability \Ai\ , which is obtained using the flux for the R function Tr, = i{hr'^)[r-'^ Rdrv"^ W - 

R'^ dj-r^ R) . As R{r) ~ U{r)/r, the usual radial wave function R is composed of incoming 
and outgoing spherical waves with amplitudes equal to the ones of the U function and the 



transmission probability jA^j is given by l-Bout/^ou 



The WKB tunnel probability 



_ 2 

Ai can therefore be directly used in Eq. to compute the value of the greybody factor. 



4.2 Results for the cross-section and radiation spectra 

Figure. [S] shows, on the left panel, the WKB emission/absorption cross-section alongside the 
exact one numerically obtained and, on the right panel, the relative error. We expect 
the WKB method to break down twice : in the low energy region because this is where the 
quantum nature of particles dominates and near the resonances, i.e lj — VJnax, because the two 
turning points are extremely close one to the other. This behaviour can be seen by looking 
at the shape of the emission/absorption cross-section. The divergence is due to the £ = 
partial wave whose transition probability tends to e~^^^^ when lo tends to zero. On the other 
hand, the cross-section at resonances, occurring when the energy of the field is equal to Vi, is 
overestimated in the WKB approach (although the spectrum is quite well recovered). The mis- 
estimation at resonances is not dramatic and leads to regular peaks in the relative error whose 
amplitude does not exceed 35%. Furthermore, it has only small consequences on the radiation 
spectrum. However, without the help of an exact calculation, the WKB approximation seems 
-as expected- hazardous to provide reliable results in the low-energy region. It is both due to 
the unavoidable behaviour of the turning points and to the huge steps for the tortoise coordinate 
numerically occuring close to the event horizon. 
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Figure 5: Left panel : Emission cross-section for a 4— dimensional Schwarzschild black hole. 
The solid line corresponds to exact numerical results, the dashed line to the first-order WKB 
approximation results and the dotted line to the improved WKB results. Right panel : Relative 
error /^ag/ag of the WKB cross-section estimation. The dashed line stands for the first-order 
WKB error and the dotted one to the improved WKB error. 



In spite of this drawback, the WKB techniques provides good results in the intermediate and 
high energy regimes. First, the high energy optical limit is well recovered, as it can be seen on 
Fig. Moreover, it was demonstrated in Ref. ^HI that if the transmission coefficient \ Ai\'' tends 
to unity when lo tends to infinity for all the partial waves, then the emission/absorption cross- 
section should tend to the optical limit when lo +oo. Our WKB transmission coefficient is in 
agreement with the above conditions and the cross-section has to tend to (Tg°'^'' = 277r7'|^/4. In 
the intermediate energy regime, the relative error admits an average value of 10% and decreases 
for higher values of the energy. 

Those results could be improved using higher order semi-classical expansions. A third-order 
WKB tunnel probability has been derived by Iyer & Will and extended to the sixth order by 
Konoplya (see Ref. Illp to derive the quasi-normal frequencies of black holes. It is optimised for 
resonant scattering near the top of the potential barrier and allows for a non- vanishing reflexion 
coefficient for particles with an energy greater than the maximum of the potential. This should 
allow to avoid the spiky structure at resonances. The transmission coefficient, successfully used 
to derive the greybody factors for gravitons emitted by 13— dimensional Schwarzschild black 
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Figure 6: Scalar radiation spectrum from a 4— dimensional Schwarzscliild black hole. The solid 
line corresponds to exact numerical results, the dashed one to the first-order WKB approxima- 
tion results, the dotted line to the improved WKB results and the dashed-dotted one to a black 
body radiation spectrum with the high energy limit of the cross-section. 



holes ^ni, takes the following form 

= i^^2M^+i/2) , ^ + 1/2 = ^P'(W) (^2^(w) j . (35) 

The greybody factors obtained with the improved WKB tunnel probability are displayed 
on Fig. [SI where it can be seen that the breakdown of the approximation at the top of the 
potential barrier is removed. However, the infrared divergence cannot be avoided with this new 
technique, as is the case for gravitons |1U) . 

The emission/absorption cross-section does not correspond to any physical quantity. The 
meaningful flux, or radiation spectrum, is displayed on Fig. |51 As it can clearly be seen, both 
WKB approximations allow for a good estimation of the high energy tail of the distribution 
whereas the low energy one is mis-estimated due to the infrared divergence. It should be 
pointed out that the WKB technique for the estimation of greybody factors is more of theoretical 
interest : it allows a first and fast estimation of the emission/ absorption cross-section in both the 
intermediate and high energy regime (which is highly relevant to understand the consequences of 
the possibly intricate metric), but does not provide any real improvement in the determination 
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of the " experimental" spectrum. 

5 Scalar fields in SAdS space-times 
5.1 Shape of the potential 

Fields in SAdS space-times have been intensively studied for quasi-normal modes estimation 
(see Ref. |11II12) and rcfcrenccs therein). Nevertheless, the general behaviour of scalar fields in 
those space-times has not been investigated and this study is carried out at the semi-classical 
order hereafter. SAdS space-times are described by a Schwarzschild-like metric with 

2M , , 

h(r) = \- — + —. (36) 

The black hole horizon radius corresponds to the positive real root of h(r) : 2M = r^f(l + 
r'j^/R'^). The potential seen by a scalar field and the tortoise coordinate are given by 



r2f^\ _ I ^ 2M I \ ( , 2M , 2 



■In 



The shape of the potential strongly depends upon rj = rn/R- If V is much smaller than one, 
the black hole horizon radius is smaller than the characteristic AdS curvature radius, while if rj 
is much greater than one. the horizon radius becomes greater than the curvature radius. The 
potential is depicted on Fig. [7| for 77 ^ 1 and is a mixing of the potentials in Schwarzschild 
space-time and in AdS space-time : when r* —00, i.e. r rn, the potential falls down to 
zero because of the event horizon whereas it diverges to -t-00 when r* — + ^g^^g, i.e. r +00, 
where rS«^q = — - — — " , ^. It also exhibits a finite barrier -because of the presence of 

the black hole- for any value of the orbital quantum number £, which leads to the appearance 

(i) . . (i) 

of a local maximum Knax and minimum V^^^^, respectively located at r*^^^^ and r*^i^. 

When 77 3> 1, the finite potential barrier is present only at high values of the angular 

momentum : as examples, the barrier appears at € ~ 10 for 7] = 1 and at £ ~ 10^ for r/ = 2. In 

the following, we will focus on the high hierarchy case (77 <C 1) since it corresponds to the most 

physically interesting one. 
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Figure 7: The potential in units of \/|A| seen by a massless scalar field in SAdS space-time as 
a function of the tortoise coordinate for th = 0.1 and for 77 ^ 1. 



5.2 Stationary states solutions 

The solution being developped in stationary states tp{r*,t) = ipr{r*)e~'^'^^ , the equation of 
motion is given by Eq. © where w is a real and positive parameter. The WKB solutions 
depends on the energy^ lu of the field and the problem has to be divided into three energy 
regimes, defined by the number of classical turning points : i) uj < V^^n (Regime I), one turning 
point located before rmL ; ii) V^^n < uj < Viiax (Regime II), three turning points, r± 
corresponding to the potential barrier induced by the black hole and rg corresponding to the 
infinite barrier at spatial infinity ; iii) Vrnlx < w (Regime III), one turning point which also 
corresponds to the infinite potential barrier. The boundary conditions, although depending 
on the considered quantity (quasi-normal frequencies, greybody factors, etc.), should read as 
follows in the most general case : 



Mr*) 



7=( 



when r* ^ — 00 
when r* rT 



(38) 



SAdS- 



Considering first the one-turning-point case, the WKB wave functions has to be considered 



^Calling energy the frequency parameter uj was meaningful in the previously studied space-times as it was 
corresponding to the energy of particles for observers located at r = in AdS Universe or at spatial infinity in 
Schwarzschild ones. Although such flat observers do not exist in SAdS background, we keep calling energy the 
frequency parameter to avoid using too many variables. 
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in two regions : i) ?' < f^/o (Region A) where the solutions are progressive waves and ii) r > 
r_/o (Region B) where the solutions are evanescent and exponentially divergent waves. With 
proper boudary conditions, the WKB wave functions are obtained by the matching procedure 
at r < r_/o described in appendix 1X1 : 



Mr*) 



^/p'{r 



{x)dx 



-i p{x)dx 



- f\ p (x)dx 



in Region A 
in Region B. 



(39) 



For the three-turning-points case, the WKB solutions are considered in four regions : i) 
r < (Region A) where the solutions are progressive waves ; ii) r_ < r < r_|_ (Region B) 
where the solutions are evanescent and exponentially divergent waves ; iii) r_|- < r < tq (Region 
C) where the solutions are also progressive waves and iiii) r > tq (Region V) where the solutions 
are only evanescent waves, taking into account the boundary conditions at ?'sAdS- -'^PPlyi'^g the 
matching procedure, the WKB solutions are found to be : 



/p(r* \ 



i p[x)dx 



-i p{x)dx 



p' {x)dx 



Air*) = < 



\/p'{r* 

Apute^^ ^gj/"" p{x)dx _ ^g- 



2^p(r* 

A2ut£llg~*f -JrJ P'{x)dx 



i p{x)dx 



in Region A 
in Region B 
in Region C 
in Region T). 



(40) 



where r = J^'^ p' {x)dx. This strange behaviour in Region B arises because we have assumed 
that the solutions can be decomposed into stationary state and because of the semi-classical 
anstaz. When considering the shape of the potential, it is clear that the states in Region 
['"max, +oo] cannot be stationary states but quasi-stationary ones, as tunneling transition from 
Region C to Region A is possible but not taken into account in a first order semi-classical 
expansion. 

The radial density probability for scalar particles can be inferred from the above stationary 
state solutions and takes the following form : 



p{r*) = i [Mrne-'^^'dt (^J (r*)e*'^*) - (r*)e'"*ft {Mr*)e''^')] = 2t^V4(^*)V'r(r*). (41) 
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Figure 8: Probability densities as functions of r* alongside the I = Q scalar potential (solid line) 
for ?; = 10 : the dashed line stands for u < V^^^, the dotted line for V^^^ < lo < V^iax and the 
dotted-dashed line for u; > vJiiax- Each probability is displayed at its energy uji, the horizontal 
line at w = corresponding to p = 0. The turning-points are shown and the energy can be 
inferred either in units of i/jAI on the left y— axis or in units of on the right. 



In the three-turning-points case, the density probability therefore reads as 



pir*) = 2uj\Aout\ I 



1 



2p(r*) 



+ cos (2 f* p{x)dx + § 
exp (-2 p'{x)dx^ 
1 + cos (2 /J^* p{x)dx + f 
exp (-2 /Jl p'{x)dx 



ip'(r*) 



in Region A 
in Region B 
in Region C 
in Region V. 



(42) 



This probability is displayed on Fig. |S1 for three different cases. As expected, at the turn- 
ing points, the probability density diverges because of the p^^{r*) factor, which is a direct 
consequence of the breakdown of the WKB approximation. Since the potential rapidly falls 
down to zero, the radial probability density rapidly tends to the coherent superposition of in- 
going and outgoing free particles, i.e. plane waves, with the same amplitude lim|.*_+_oo p(^*) = 
^ [1 + cos (2ijjr* + f + '/'(''fni))] 7 where (/? is an additional phase depending on the initial posi- 
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tion chosen for integration. The minimum and maximum energy cases can be divided into two 
parts : before the turning point, the probabihty density is described by interferences between 
ingoing and outgoing modes, whereas after the turning point the probabihty density is a de- 
creasing exponential which faUs down to zero at r* = '"g^ds- I^^r the intermediate energy case, 
the space is separated into four regions : for r < r- and r+ < r < ro, the probabihty density 
takes the form of interferences between ingoing and outgoing modes whereas it takes the form 
of a decreasing exponential in the r_ < r < and r > tq regions and falls down to zero at 

'''* — ^SAdS- 

Some conclusions about the black hole evaporation process with a negative cosmological 
constant can be inferred in this framework. As in SdS space-time, the evaluation of the trans- 
mission coefficient critically depends upon the position of the observer. If the observer is located 
at spatial infinity, the transmission coefficient falls down to zero because of the total reflexion 
at rg. However, this can be avoided by considering an observer near the local minimum, which 
is always possible in the case of a high hierarchy. Particles with an energy lying in Regime I 
are then totally reflected at r_ and the observer should not see any particle with an energy 
lower than Vjfjjjj coming from the black hole. In addition, for quanta with lo > V^^l^, only modes 
with angular momentum lower than to, so that V^™ < a; < V^^^ , can reach the observer. 
For Lu < Vmax, the transmission coefficient is given by e"^"^ while it becomes equal to one for 
particles with u > Vmax- The qualitative features of the evaporation process of SAdS black 
holes can therefore be captured. It should, of course, be pointed out that a realistic study 
of this phenomenon has to be carried out dynamically, involving a full resolution of the KG 
equation while keeping the time dependence. 

5.3 Energy levels and band widths 

Because of the shape of the SAdS potential, the density of states for scalar fields is smooth 
and continuous with some resonances at specific energies, lying in the [VJnin , VJnax] range, with 
bandwidths given by the tunnel probabilities. This should lead to stationary states when 
the potential barrier becomes either infinitely large or infinitely high. The aforementioned 
conditions can be obtained for high angular momentum quantum numbers as well as for a 
vanishing rj ratio. Working at the semi-classical order, the energy levels are computed by 
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n = Q 


n = 1 


n = 2 


n = i 


^ = 


2.49-^0.24 


4.51 - iO.44 


6.52-^0.63 


8.53-^0.83 


I = 1 


3.82-i0.02 


5.82 - i0.04 


7.82- i0.05 


9.82-i0.07 


1 = 2 


4.86- il.6 10-^ 


6.86-^2.3 10-^ 


8.86-^3.0 10-'^ 


10.85 -i3. 7 10-^ 


^ = 3 


5.88- i9.0 10-^ 


7.88-^1.2 10"* 


9.88-^1.5 10-4 


11.87 -il. 8 10-"* 


^ = 4 


6.88- i4.6 10-« 


8.88-^6.0 10-« 


10.88 -i7.3 10-" 


12.88 -i8. 7 10"*^ 


i = b 


7.89-^2.3 10- 


9.89 -i2. 9 10-" 


11.88-^3.5 10-" 


13.88 -i4.1 10"" 



Table 1: Complex energy levels (vn,i x R) for 77 = 5 x 10 ^ as functions of the angular momen- 
tum quantum number £ and the integer n of the Bohr-Sommerfeld quantization rule. 



solving numerically the following algebraic equation : 



2 / — dr - (2n + 1) TT = 0, (43) 

where ^ stands for the unknown quantity. Unlike the AdS case, the usual radial coordinate r 
is prefered over the tortoise one. The numerical investigations show that, at the semi-classical 
order, the product w„.£ x R only depends on the ratio 77, (n, €) being fixed. Consequently, the 
energy levels can be written in a very convenient way : ujn/ x i? w f{n,£) x 5(77), where / 
and g are two functions which respectively depend on (77,, £) and 77. The bandwidth T of the 
resonances is given by the product of the timnel probability \Ai\^ with the number of collisions 
at the turning point r+ : Tnj = u!n,e x \A(\^ , with \A(\^ given by 




lAil"^ = < V / (44) 

(^1 + cxp (^-2V(w) (2^(w)) '^'^^ • 

In view of the results obtained for the greybody factors, the second WKB tunnel probability 
should be preferred as it provides more accurate results. The numerical investigations lead to 
assume that this bandwidth can also be separated into a (77, £) dependence and a 7/ dependence, 
as in the case of the resonant energy position. 

The semi-classical characteristics of resonances numerically obtained are presented in Ta- 
bles ^ to 01 for four different values of the ratio 77 = {0.0005,0.001,0.01,0.1} corresponding to 
cases where the finite potential barrier appears at any multipolar order. The tables give the 
complex energy i = uin^e — i^n.e in units of i?-^. As the KG equation is solved with the 
ansatz = e~*''*Y7"(n)i?(r), the imaginary part of the complex energy has to be negative. 
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n = Q 


n = l 


n = 2 


n = 3 


^ = 


2.51 - iO.24, 


4.54 - iO.44 


6.55-^0.64 


8.56-^0.83 


1 = 1 


3.83-^0.027 


5.82 - i0.042 


7.82- i0.057 


9.81 - t0.071 


1 = 2 


4.86- il.6 10-^ 


6.86-^2.3 10^^ 


8.86-^3.0 10-'^ 


10.85-^3.7 10-^ 


i = i 


5.88-^9.0 10-^ 


7.87 -il. 2 lO-'' 


9.87 -il. 5 10-4 


11.87 -il. 8 lO-'' 


i = A 


6.88- i4.6 10-« 


8.88-^6.0 10-" 


10.88 -i7.4 10-'^ 


12.88-^8.7 10-« 


£ = 5 


7.89- i2.3 10- ^ 


9.89-^2.9 10-^ 


11.89 -i3. 5 10-^ 


13.88-^4.1 10-^ 



Table 2: Complex energy levels {vn,£ x R) foi' ?/ = 10 ^ as functions of the angular momentum 
quantum number £ and the integer n of the Bohr-Sommerfeld quantization rule. 





n = 


n = 1 


n = 2 


n = 3 


£ = 


2.61-i0.26 


4.64 - iO.47 


6.63- i0.70 


8.62-i0.97 


£ = 1 


3.82-i0.028 


5.79- i0.044 


7.77- iO.060 


9.74- i0.079 


£ = 2 


4.85- il. 7 10-^ 


6.83-^2.4 10-^ 


8.81 - i3.2 10-'' 


10.78 -i4.1 10"^ 


£ = 3 


5.87- i9.2 10-^ 


7.85-^1.2 10"* 


9.83-^1.6 10--* 


11.80 -i2.0 10-* 


£ = 4 


6.88- i4.8 10-« 


8.86-^6.2 lO-'^ 


10.84 -i7.7 10-** 


12.82 -i9. 3 10-« 


£ = 5 


7.88-^2.4 10-"^ 


9.87-?:3.0 10-' 


11.85-^3.7 10-^ 


13.83 -i4.4 10-"^ 



Table 3: Complex energy levels {i'n,i x R) for 77 = 10 ^ as functions of the angular momentum 
quantum number £ and the integer n of the Bohr-Sommerfeld quantization rule. 





n = 


?i = 1 


n = 2 


n = 3 


£^0 


2.65-i0.75 


/ 


/ 


1 


£ = 1 


3.6- i0.12 


5.4-^0.95 


/ 


1 


£ = 2 


4.74- i8.6 10-^ 


6.52- i0.055 


8.22- iO.459 


9.81 - i3.42 


£ = 3 


5.77-^5.5 10-* 


7.58 -i2. 8 10-^ 


9.34- i0.019 


11.03 -iO. 15 


£ = 4 


6.79-^3.4 10-^ 


8.62 -il. 5 10-* 


10.41-^8.2 10-* 


12.15 -i5. 5 10-^ 


£ = b 


7.80- i2.0 10-** 


9.65-^8.3 10-*^ 


11.46 -i3.9 10-* 


13.23 -i2. 2 10-* 



Table 4: Complex energy levels (vn^i x i?) for 77 = 10 * as functions of the angular momentum 
quantum number £ and the integer n of the Bohr-Sommerfeld quantization rule. 
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The positions of resonances tend to the pure AdS energy spectrum when rj tends to zero for 
aU angular momenta except when £ = 0. This is not surprising since, for SAdS space-times, an 
infinite steepness of the waU never arises, even for arbitrarily small black holes : the transition 
from AdS to SAdS space-times is singular. 

The bandwidth of the resonances behaves as expected : on the one hand, either a decrease of 
77 or an increase of £ leads to a decrease of the bandwidth since the well-potential is deeper and, 
on the other hand, the bandwidth becomes larger for higher values of n because of the increasing 
tunnel probability. Finally, when 77 ^ 0, the quasi-stationary states tend to the stationary ones 
of an AdS Universe except for a vanishing angular momentum quantum number. The specific 
behaviour of the monopole partial waves underlines the singular behaviour of the SAdS scalar 
potential. However, it cannot be excluded that this arises from the use of semi-classical methods 
and is only an artifact of this approximated scheme. 

6 Conclusions and perspectives 

It was shown in this article that semi-classical methods for the resolution of Schrodinger-like 
equations can be successfully applied to scalar fields propagating in a static, spherically sym- 
metric, curved background. The radial part of the KG equation has to be rewritten in a more 
convenient in a systematic scheme : the tortoise coordinate is used and the unknown func- 
tion is redefined independently of the metric. The propagation of fields is then described by a 
1— dimensional Schrodinger equation whose resolution can be performed with a semi-classical 
ansatz. However, the classical problem underlying the quantum one is based on a slightly- 
modified-from-relativistic-case Hamilton-Jacobi equation where what is usually interpreted as 
the potential is here seen as the square of the potential. The semi-classical method was then 
used to study the behaviour of scalar fields in AdS and Sehwarzschild backgrounds and it was 
shown that the energy levels in an AdS universe and the greybody factors entering the Hawking 
radiation law are both recovered at the semi-classical order with a very good accuracy. It should 
also be emphasized that the WKB estimation of the greybody factors provides a powerful tool 
for theoretical investigations, i.e for testing the influence of a given space-time curvature on the 
emission/absorption cross-sections. Finally, we have performed a semi-classical study of scalar 
fields propagating in SAdS space-times. The potential seen by those fields imposes resonances 
in the density of states, whose positions in energy and bandwidths can be evaluated using the 
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Bohr-Sommcrfeld quantization rule for the former quantity and the WKB tunnel probability 
for the latter one. 

In principle, all the above results could be improved : first, the Taylor expansion of the 
ansatz used in the KG equation could be preformed up to higher orders. Then, as it has 
already been mentioned by Langer considering the radial part of the Schrodinger equation 
as a one dimensional equation remains a naive approach as the centrifugal potential exhibits a 
singular behaviour when r ^ 0. As discussed in ^^Ij the phase of the WKB wave function for 
radial motion docs not match the one of the exact wave function in regions where the potential 
can be neglected when compared to the centrifugal one. Nevertheless, the matching between 
both solutions can still be performed by replacing £{£ + 1) by {£ + 1/2)^ but, in any case, a 
precise study of the matching between the semi-classical wave functions and the exact ones, in 
regions where such exact solutions can be found, have to be carried out before applying the 
Langer trick. 

Finally, the semi-classical methods provide a powerful tool to study the dynamics of scalar 
fields in curved space-times as, in most cases, exact solutions cannot be derived. As an example, 
the highly non-trivial problem of the Hawking evaporation process for SAdS black holes can 
be greatly simplified in this way and is currently the only available treatment. Another field of 
application is related to primordial cosmology where the WKB approach, first used by Martin 
& Schwarz jSj, has been extended to the second order by the authors of Rcf. [HI so as to 
derive the primordial fluctuations power spectrum. However, it has been recently proven by 
the authors of Ref. ^2 that the accuracy of the WKB approximation is intrinsically limited 
because it involves divergent series. This result, demonstrated in the case of cosmological 
particle production during the inflationary era, enlightens the cautious needed when dealing 
with WKB techniques. 
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A Appendix IXt Matching procedure 



The matching procedure involves the equation + p'^{0''P{0 = with two turning point 
at ^1 and ^2- The procedure consists in restricting the momentum at the linear order p^(^) ~ 
— A(^— ^i), with A > 0, and to match the new solutions with the WKB ones on the left and on the 
right of the turning point [51l51ll5|. With this linearized momentum, the solution of the previous 
equation are expressed in terms of Airy functions (Ai, Bi) : ip{^) = -BiAi ^A^ — ^i)^ + 
B2^i (^A3(^ — ^i)^ . As ^ is close to ^i, this solution can be expanded in terms of exponential: 



m = { 



X 12 



2v^ 



^ - CiP^ [(^2 - ^Bi) exp (zf VA - 



+ i'. 



(S2+^Sl)exp(-^|^/A|e-aP -«f 



A~T3 



+B2exp(|VA|^-a 



f exp (-I^AI^-a 



in region A 



(45) 



in region B. 



To perform the matching between Airy functions solutions and the WKB solutions, we calculate 
the WKB solutions in the two region using the linear approximation for the momentum : 



m = { 



A--ie-ei 



Aout cxp (-i^VXl^ ~ 6 



A„icxp (i|\/A - Ci| 



in region A 



(46) 



^"^1^-61 ^^e«;aexp(-|VA|^-ai') in region R 



Matching i?i/2 to Aout/in and to Aeua/diuj it is possible to match the two solutions in regions 
A and B : Aeva = e^'^ Aout and Ai„ = e~^^ Aout- 

The same procedure is applied at the second turning point ^2 just by introducing ^' — — ^, 
at this point the WKB solutions take the form : 



V'(S') 



A- 



le' - ?2r^ A,uae-^ exp (|VA - ^^p) in region B 



(47) 



A~^ w - e^r^ Bout exp Uivxw - e^p 



in region C 



with T = p'iO'^^ and gives: Bout = e Aeva- From this equation and the previous one, 
we finally obtain the amplitude of the outgoing waves at +00 in function of the amplitude of 
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the outgoing waves at — oo: Bout = e '^Aout- 

B Appendix El Resolution of the AdS-KG equation 

To solve |5] the differential Eq. (|f 9|l . a change of variable as well as a change of the unknown 
function is performed : z = cos (2p) ; R{p) = s'm^ (p) coa^ {p)P{p). Under this change, the 
differential equation is of the form 15" {l-z'^)^ + {b-a-{a+b+2)z)^+n(n+a+b+l)P(z) = 
where the parameters (/?, a, b, n) must satisfy the four constraints : 

a + b + 2^e+l + P; b-a = f3-2~e; /?(/?- 3) = 0; n{n + a + b + I) = I {ly^ ^ {£ + P)^) . 

(48) 

The solution are the Jacobi's polynomials P^'^\z), imposing that n is an integer, a > —1 
and b > —1 ^|. Once the four constraints for the parameters {(3,a,b,n) solved, the solution 
of the radial part of the KG equation is expressed in function of the Jacobi's polynomials : 
i?(p) ~ sin^ (p) cos'^ (p)^(^,!-^-3)%^^ ('^^^ where we choose (3 — 3 in order to respect the 
6 > — 1 constraint. 

C Appendix O Determination of the turning points 

In the £ Q case, the determination of r^. is more complicated. Introducing y = tan {r* /R), the 
solutions of the equation ui = V{r) are given by the algebraic equation 2y'^+(2 + £{£ + 1) — upR^^ 
£{£ + 1) = 0. In order to have two different, positive, real roots, three constraints have to be 
satisfied: 

2 + £{£+!)- uj'^R^ <Q 

{2 + £{£+l)-uj^R^f -U{£+l)>Q (49) 
8^(^+1) > 

where the third one is always satisfied. Solving the first constraint, the zero-point energy feature 
of the energy spectrum is recovered w > \/2 + £{£ + 1). This result is in good agreement with 
the exact energy spectrum as y^2 + £{£ + 1) < £ + 3. Under such constraints, only two positive, 
real roots exist which are the two turning points: 

y± = ]^\juj^R'^ -2-£{£+l)±\l{2 + £{£+l)- uj^K^f - 8£{£ + 1). (50) 
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